Superconductivity in carbon nanotube ropes: Ginzburg-Landau approach and the role 

of quantum phase shps 
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We derive and analyze the low-energy theory of superconductivity in carbon nanotube ropes. A 
rope is modelled as an array of ballistic metallic nanotubes, taking into account phonon-mediated 
plus Coulomb interactions, and Josephson coupling between adjacent tubes. We construct the 
Ginzburg-Landau action including quantum fluctuations. Quantum phase slips are shown to cause 
a depression of the critical temperature Tc below the mean-field value, and a temperature-dependent 
resistance below Tc. 

PACS numbers: 73.63.Fg, 72.10.-d, 74.25. Kc 



Over the past decade, the unique mechanical, electri- 
cal, and optical properties of carbon nanotubes, including 
the potential for useful technological applications, have 
created a lot of excitement. While many of these prop- 
erties are well understood by now, the experimental ob- 
servation of intrinsic M 0, 0] and anomalously strong 
proximity-induced 0, [3 superconductivity continues to 
pose open questions to theoretical understanding. In this 
paper we analyze ID superconductivity found in ropes 
of carbon nanotubes P, Q , starting from a microscopic 
model of the rope in terms of an array of individual clean 
single- wall nanotubes (SWNTs), with attractive phonon- 
mediated on-tube interactions and intcr-tube nearest- 
neighbor Josephson couplings. Then the Luttingcr in- 
teraction parameter gc and the Josephson coupling A 
describing transfer of Cooper pairs are crucial micro- 
scopic parameters, which^fortunately, can be estimated 
rather accurately 0, Q, I3- The coupled-chain prob- 
lem corresponding to superconducting nanotube ropes, 
where typically less than hundred metallic SWNTs are 
present P, 0], does neither permit classical Ginzburg- 
Landau (GL) theory nor a standard self-consistent BCS 
approach, in contrast to the situation encountered in, 
e.g., quasi-lD organic superconductors |^. At this time, 
nanotube ropes represent wires with the smallest number 
of transverse channels showing intrinsic superconductiv- 
ity, even when compared to the amorphous MoGe wires of 
diameter « 10 nm studied in Ref. [l^ , where still several 
thousand channels are available. Based on a microscopic 
derivation of the quantum GL action, we show that quan- 
tum phase shps (QPS's) [ill [iS 13 are crucial for an 
understanding of experimental results 0,13- First, they 
cause a depression of the transition temperature Tc below 
the mean-field critical temperature T^. Furthermore, for 
T < Tc, a finite resistance R{T) due to QPS's appears, 
which exhibits approximate power-law scaling. Below we 
determine the full temperature dependence of R(T < Tc) 
for arbitrary rope length. 

We consider a rope consisting of N metallic SWNTs, 
where disorder is assumed to be negligible '25^. The va- 
lidity of modelling the rope as an array of ballistic ID 
quantum wires has recently been discussed in Ref. 
Since the K point degeneracy is inessential here, an in- 



dividual SWNT can be described as a spin-1/2 Lut- 
tingcr liquid, where the combined effects of Coulomb and 
phonon-mediated interactions lead to an interaction pa- 
rameter gc 3|, where gc = 1 refers to the nonintcracting 
case and 5c > 1 {9c < 1) signals effectively attractive (re- 
pulsive) interactions. We mention in passing that lattice 
commensurabilities and electron-electron backscattering 
can be neglected in the intrinsically doped SWNTs en- 
countered in practice |l5l llSj. In a thick rope, Coulomb 
interactions are expected to be largely screened off, and 
5c > 1 due to breathing-phonon exchange The only 
weak screening in the thinnest ropes studied in Refs. P,y| 
is probably linked to the absence of superconductivity in 
these samples. To probe superconductivity in ultrathin 
ropes, it is necessary to externally screen Coulomb inter- 
actions. In principle, three different inter-tube coupling 
mechanisms should now be taken into account, namely 
(i) direct Coulomb interactions, (ii) Josephson coupling 
(Cooper pair hopping), and (iii) single-electron hopping. 
The last process is strongly suppressed due to the gen- 
erally different chirality of adjacent tubes and, in 
addition, for gc > 1, inter-SWNT Coulomb interactions 
can be neglected Therefore the most relevant mech- 
anism is Josephson coupling between adjacent SWNTs. 

In the (idealized) rope crystal, (metallic) SWNTs of ra- 
dius R are arranged on a trigonal lattice with lattice con- 
stant a — 2R + b, where b — 0.34 nm JJj. Choosing the 
X-axis parallel to the rope, and numbering the SWNTs by 
j — 1, . . . , N, with center at fj = {yj, zj) = niEi + 71202, 
where ai = a(l, 0) and 02 =a(l/2, a/3/2) span the trigo- 
nal lattice, allowed indices (tt-i, 712) corresponding to j fol- 
low from the condition \rj\ < i?ropc- A given rope radius 
i?i.opc then fixes the number of tubes N. The Josephson 
coupling matrix Ay is nonzero only for nearest-neighbor 
pairs (?,j), where A^ = A. Only singlet pairing of elec- 
trons on the same tube is important leading to the 
on-tube order parameter Oj{x,T) = J2r<T '^''{'ro-j4'-r,-<Tj, 
where ipraj{x,T) is the electron operator for a right- or 
left-moving electron (r = ±) with spin cr = ± on the 
jth SWNT, and < r < 1/T is ima gina ry time (we put 
S = fcs = 1). In bosonized language 
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Oj — (Trao) ""^ exp 
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where gq = 0.246 nm is the SWNT lattice spacing, 
'^c/s,j(2;, t) denotes the charge/spin boson field on the 
jth SWNT, and ^c/sj is the dual field to ipc/s,j- The 
Euclidean action is 

N 

S^J2 "^LLi^cj, ^s,,] - J2 J dxdTO]0„ (2) 

J = l jk 

where the on-tube fluctuations are governed by an effec- 
tive Luttinger liquid action [sl ITHlll^. 



'5'lL [Sc, Vs] 



dxdr 



+ 



with Vc/s — vp/gc/s for Fermi velocity vf — & x 
10^ m/sec, and gs ^ 1 due to spin SU{2) invariance. 
Note that the model (0) and our results below apply 
beyond the specific system under study here, see also 
Refs. Hill. 

Next we employ a Hubbard- Stratonovich transforma- 
tion in order to decouple the Josephson term in Eq. |(2Jl, 
using the complex- valued order parameter field Aj(a;, r). 
This allows to write the partition function as 

j VAcxp (^-Y,So[A,]- J dxdT^A*Z?,feAfcj , 

(3) 

where A* is the complex conjugate field. The N x N 
matrix D denotes the positive definite part of A^^, since 
order parameter modes corresponding to negative eigen- 
values of A can never become critical. Integration over 
the on-tube boson field fiuctuations then leads to 



5o[A] = - In / Vej)ipse 



-SLLlfcVs]-/ da;(ir(A'C'+AOt) 



(4) 

The expectation value of Eq. ^ can be computed as 
{Oj) = ^k-^jk{^k), where (A^) is averaged using the 
action corresponding to Eq. 

Assuming a spin gap, taking T = 0, and allowing 
only static homogeneous configurations Aj(x,T) — Aq, 
Eq. Q can be evaluated explicitly 19]. Saddle-point 
analysis then yields a relation between the mean-field 
critical t emp erature and the T = superconducting gap, 
see Ref. [l^. For general order parameter Aj{x, r) or ar- 
bitrary temperature, however, integration over the Lut- 
tinger phase fields in Eq. Q is impossible. To make 
pro gres s, it is instructive to construct the GL action 
|l3l |2(]| |. where it is crucial to include quantum fluctu- 
ations. A systematic approach proceeds via cumulant 
expansion of Eq. Q up to quartic order in the expansion 
parameter \A\/2ttT j2fl]. We stress that this expansion 
is carried out for the single-chain problem, and is not 
restricted to N ^ 1 . Assuming slowly varying conflgura- 
tions Aj(a;,r), gradient expansion yields the Lagrangian 



in Eqs. © and as 

L = J2^,,A*A, + J2{m-'-A]\A,\' (5) 

+ B|A,|4 + C[|a.A,f + c;2|a.A,f]), 
with the Mooij-Schon plasma velocity [23|, Cs = 



Vc\JC/b, and V,k - Ea(W^a W^r')01«>(a|^>, where 
Wa denote the eigenvalues of A in descending order (a = 
1, • • • tN), with eigenvectors |a), and the a-summation 
extends over Wa > only. Furthermore, the positive 
coefficients A, B, C are with 7 = Vc/vg given by 



A 



B = ^:^^{naoT/v.f^^'+'^^-'B, 



c 



327r4w3 



(6) 



Putting gs = 1, dimensionless ^c-dependent numbers are 
defined as 



C 



dz- 



w 



(7) 



fc{z)fAzy 

where we use the notations z = (w,u), 

fc{z) = I sinh(w + fs{z) = \ sinh(7W + iu)\, 

and J dz = du dw. Here A (D) is given by Eq. {T)) 
with ^ 1 (w^ u^), and 



B = 



dzidz2dz3 

.fc(z2)fc{zi3) 
fs{zi)fs(z23) 



fs{z2)fs{zi'i)fs{Z'i)fs{zi 



4 fc{zi)fc{z23) 

fs{z2)fs{zi3) fc{z3)fc{zi2) 

-I- (1 ^ 2) + (1 ^ 3) 



with Zij = {wi — Wj,Ui — Uj). The quantity B is eval- 
uated using the Monte Carlo method, see also Ref. 
For gc = I, we first numerically reproduced the exact 
result 13 = Stt^C with C = 77rC(3)/4 Numerical 
values can then be obtained for arbitrary gc- ^From 
Eq. (O, previous GL results for the infinite 2D array 
of coupled ID chains are recovered |^ . In that case, 
the Vij term in Eq. jSJ) leads to transverse gradients be- 
cause a corresponds to transverse momentum kj_, with 
W-^ - W^^ cx P^. Note that Eq. © additionally in- 
cludes quantum fiuctuations and allows to describe the 
case of arbitrary TV. 

^From Eq. (|SJ), we obtain the mean-field critical tem- 
perature 
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(8) 



3 



0.2 r 



0.15- 



H 

0.1- 

o 

<l 



0.05- 



phase slips ■ Following arguments similar to the ones 
of Ref. we find that only QPS's play a role. Nu- 
merical evaluation of Eq. (fTn|l shows that well below T°, 
transverse fluctuations are heavily suppressed [2^, and 
therefore QPS's can be described using the action 



(12) 



where /i — /ij is a dimensionless rigidity. For not too 
low temperature, and neglecting transverse fluctuations, 
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^l{T) = aoN \l - {T/T°Ysa-i)/g. 



(13) 



FIG. 1: Temperature dependence of Ao/2ttT versus T/T° 
for N ^31 (open) and iV = 253 (filled circles). 



Assuming sufficiently thick ropes such that Coulomb in- 
teractions can be neglected, in concrete estimates we 
shall put 5c = 1-3 jSI], with Josephson coupHng X/vp = 
0.02 Numerical evaluation yields 

i~ 30.72, i?~ 293.1, 12.12, Z) ~ 7.78. (9) 

Equation © then predicts, e.g., T° = 2.3 K for TV = 
31, which is slightly above reported experimental values 
fn what follows, we focus on temperatures below 
T°. Writing Aj — |Aj| exp{i0j(a;, r)}, the amplitudes 
|Aj| are finite, with a gap for fluctuations around the 
mean-field value. At not too low temperatures, they are 
found from the saddle point equation 



where ao — 47rA((7I))^/^/i?, resulting in ao — 12.7 for 
gc = 1-3. Remarkably, at T = 0, Eq. (|13|l coincides up to 
a prefactor of order one with the rigidity p, obtained from 
standard mean-field relations _2Q|, fl — tt'^UsR'^ /2m*Cs — 
aoN, where Ug is the density of condensed electrons. At 
T — Q, this implies ao ~ vp/cg- We conclude that the 
GL prediction (|13() for //(T) is robust and useful even 
outside its validity regime. 

QPS's are topological vortex-like excitations of the 
superconducting phase field 0(a;,r). For rope length 
i — > oo and thermal length Lt = Cs/ttT —> oo, a 
QPS with core at {xi,Ti) and winding number Vi = ±1 
(higher winding numbers are irrelevant) is (f>{x, t) = 
Vi tan~^[(a: — Xi) / Cs{t — r^)] p^ . where the finite L, Lt 
solution follows by conformal transformation. This form 
solves the equation of motion for Eq. (|12|l with a singu- 
larity at the core, where superconducting order is locally 
destroyed. The local loss of condensation energy density 
Ec (this may also contain other energy costs [ij) leads 



^ Vij\Aj\ + (T4^f ^ - A)|A,| + 2B|A,p = 0, (10) to the core action Sc = k^Ec/cs, with core radius k as 



whose numerical solution (via a Newton-Raphson root 
finding scheme) yields the transverse order parameter 
profile and Aq = ^ |Aj|/iV. Typical results for A^j^ixT 
are shown in Fig. ^ which demonstrates that GL theory 
quantitatively holds down to T « 2^c/2- In our discus- 
sion below, it is useful even down to T = 0. Fixing the 
amplitudes |Aj| at their mean-field values, the result- 
ing Lagrangian governing the massless phase fluctuations 
(Goldstone modes) is 

^ = Y.^Vsid^'^i)' ^c-\d^<t>,f\ (11) 



^2y,,|A,||A,|cos( 



i>j 



with dimensionless quantities = 27rC|Aj p/cs. Elec- 
tromagnetic potentials can then be coupled in by stan- 
dard Peierls substitution rules e.g. allowing to de- 
scribe the Meissner effect. Furthermore, dissipative ef- 
fects can be included following Ref. [T^ . 

In the ID situation encountered here, superconductiv- 
ity can be destroyed by thermally activated or quantum 



variational parameter. The optimal value of the core ra- 
dius is K = (csiJ-f^EcY^^ , where Sc — fJ,/'2, and k now 
serves as UV cutoff length of the field theory. To lead- 
ing order in k/ L, k/ Lt, the hydrodynamic action l|12|) 
of a vortex is Sei = /iln[min(L, Lt)/2k] -I- S'{L/Lt), 
where L/Lt measures the anisotropy of this finite-size 
2D Kosterlitz-Thouless problem. In particular, we find 
S' {Lt ~ L) ~ 0.11/i, while in the opposite limits, 
S" ~ 2/i/7r. Below S' is taken into account as renor- 
malization of S'c. 

The next step is to analyze a QPS gas, where text- 
book analysis |23| leads to the picture of an interacting 
Coulomb gas of charges Vi = ±1, fugacity y = e~"^'=, 
and total charge zero. For fjL > ^i* k, 2, QPS's are con- 
fined into neutral pairs, quasi-long range superconductiv- 
ity is present, but QPS's cause a finite resistance below 
|11|. For fl ^ fi* , QPS proliferation leads to a Kosterlitz- 
Thouless transition to the normal metallic state. (Of 
course, here "normal" does not imply Fermi-liquid be- 
havior.) The transition temperature Tc is therefore not 
T° but follows from the condition /i(Tc) = fj,* . Equation 
((T5)l then yields 



T,/T° = [1 - f,*/aoNr^^'^-'^ 



(14) 
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This Tc depression is normally rather weak, e.g. for 
= 31, we obtain T^/T^ = 0.97, but for small N, the 
effect can be large. Furthermore, other mechanisms not 
included in our model could act to effectively reduce ao 
and hence Tc, e.g. disorder and heating effects or 
the electromagnetic environment [T^ . 

The temperature dependence of the linear resistance 
R{T) = V/I for T < Tc can be obtained by comput- 
ing the voltage drop V for applied current /. Expanding 
the vortex partition function up to order and extract- 
ing the imaginary part of the free energy F(I) using the 
Langer approach [2J, r(±/) = -2ImF(±/) can be in- 
terpreted as the rate for a phase slip by ±2tt ^J. The 
average change in phase is then (0) = 27r[r(J) — r(— /)], 
which from the Josephson relation implies a voltage drop 
V = n[T{I) - r(-/)]/e. Using e = Trhl/e, the rate r(e) 
follows for L, Lt S> k but arbitrary L/Lx in the form 



J-L/2 



^^^iet-fj.[G{t+x/c^)+G{t-x/c,)] 



where G{t) = In [(Lt/k) sinh(7rr|t|)] + i(7r/2)sgn(t). To 
evaluate the rate r(e) for arbitrary L/Lt, we replace 
the boundaries for the x-integral by a soft exponential 
cutoff, switch to integration variables t' = t — x/cg and 
t" ~ t -\- xj and use the auxiliary relation 



\t" -t'\IL _ _£f_ 



ds- 



-is{t" -t') 

+ {cslLf 



It is now straightforward to carry out the t' ,t" time in- 
tegrations, and some algebra yields the linear resistance 



Rg \T{fi/2 + l/2) 



(15) 



du- 



2/tt 



r{p/2 + iuLT/2L) 



r(M/2) 



in units of the resistance quantum Rq = Trh/2e'^. Equa- 
tion H15(l leads to good agreement with experimental data 
0,13; a detailed comparison will be given in Ref. [23|- For 
L/Lt 3> 1, the It- integral approaches unity, and hence 
R cx T'^f^-^^ while for L/Lt ^ 1, dimensional scaling ar- 
guments give R (X T^^^^. The exponents are determined 
by the temperature-dependent stiffness H13|l . While both 
power-law behaviors have been reported in Ref. |ll|. 
Eq. (|15|l describes the full crossover for arbitrary L/Lt- 
In Refs. 0,01, typical lengths were L « l/xm, which in- 
deed puts one into the crossover regime Lt ~ L. 



To conclude, we have studied superconductivity in car- 
bon nanotube ropes, starting from a model of ballis- 
tic SWNTs with attractive intra-tube interactions and 
inter-tube Josephson coupling. We have constructed 
the Ginzburg-Landau theory including quantum fluctua- 
tions. This allows for detailed predictions about the crit- 
ical temperature Tc and the QPS-induced resistance be- 
low Tc- If repulsive Coulomb interactions can be screened 
off efficiently, our theory suggests that superconductivity 
may survive down to only a few transverse channels in 
clean nanotube ropes. — We acknowledge useful discus- 
sions with A. Altland, H. Bouchiat, F. Essler, and A. 
Tsvelik. This work has been supported by the EU net- 
work DIENOW and by the SFB-TR 12 of the DFG. 
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